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Pauli spin blockade (PSB) is a significant physical effect in double quantum dot (DQD) systems. In 
this paper, we start from the fundamental quantum model of the DQD with the electron-electron in¬ 
teraction being considered, and then systematically study the PSB effect in DQD by using a recently 
developed non-perturbative method, the hierarchical equations of motion (HEOM) approach. The 
physical picture of the PSB is elucidated explicitly and the gate voltage manipulation is described 
minutely, which are both qualitatively consistent with the experimental measurements. When dot- 
dot exchange interaction is involved, the PSB effect may be lifted by the strong antiferromagnetic 
exchange coupling. 


I. INTRODUCTION 

Pauli spin blockade (PSB) is an important physical 
effect appears in double quantum dot (DQD) systems, 
which was discovered experimentally in vertically cou¬ 
pled GaAs/AlGaAs DQD as early as 2002 [^. The basic 
picture is that the hopping of electrons between two dots 
will be influenced by their spin configuration if the total 
excess electron number of the system is 2 with occupation 
state {Ni,N 2 ) = (2,0), (1,1) or (0,2), as a consequence, 
the current - voltage(/—U) curve will show a rectification 
behaviour. Obviously, the PSB is caused by the universal 
Pauli’s exclusion principle. It receives extensive studies 
in various quantum dot systems with different structures 
from vertical to lateral dots Q and from double to three 
dots 3, as well as in dots with different semiconductor 
materials from GaAs/AlGaAs to Si Recently, the 
PSB has been used to fabricate and readout the singlet- 
triplet spin qubit, which willpromote the development 
of the quantum information [^. 

Some important characters in the PSB regime have 
been investigated by virous theoretical groups. Those 
include: 1) the correlation between the PSB effect and 
occupation of the two-electron triplet state 2) the 
dynamical nuclear spin polarization by hyperfine interac¬ 
tion 0 ; ^ the nonthermal broadening effect of tunneling 
current [Sj; 4) the leakage-current line shapes from in¬ 
elastic cotunneling 5) the spin-flip phonon-mediated 
charge relaxation in double quantum dots 0; and 6) 
the PSB and the ultrasmall magnetic field effect in or¬ 
ganic magnetoresistance 0. The Pauli master equation 
(PME) with second order Fermis golden rule is the main 
approach in above works to archive the transition rates. 
Other approaches (such as nonequilibrium Greens func¬ 
tions) are not yet so popular in literatures [l^ - [l^ . 

We would like to comment that the PME is not ac¬ 
curate enough for the PSB theory. Firstly, the DQD 
is a typical quantum open system with infinity degree of 
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freedoms of the total density matrix, while the PME only 
concerns the diagonal terms of the reduced density ma¬ 
trix and treats the dot-electrode couplings by low-order 
perturbation schemes; and secondly, the DQD is also a 
typical strongly correlated system with infinity degree of 
freedoms of the electron-electron (e-e) interactions, while 
the PME either neglects this important interaction or 
treats it in the single electron level. 


Obviously, for the theoretical study on such fundamen¬ 
tal physics processes as the PSB, a non-perturbative ap¬ 
proach is highly required to deal with the basic quantum 
model involving the e-e interactions. The hierarchical 
equations of motion (HEOM) approach we newly devel¬ 
oped can meet this requirement, which nonperturbatively 
resolves the combined effects of dot-electrode dissipation, 
e-e interactions, and non-Markovian memory [15l - l22j |. In 
this paper, we start from the Anderson multiple impurity 
model to describe the DQD, fully considering the e-e in¬ 
teraction and the dot-electrode couplings. By using the 
HEOM approach, we deal with this quantum model non¬ 
perturbatively to accurately obtain some observations, 
such as the spectral function, occupancy of electron spin 
and current, etc. Our theory not only can reveal the 
physical picture of the PSB clearly but also can eluci¬ 
date its dependence on various parameters such as the 
gate voltages, dot-dot coupling, and exchange correla¬ 
tion between spins in different dots. Besides, the exter¬ 
nal field manipulation is also convenient to be involved 
in our theory. 


The paper is organized as follows. In Sec.H we briefly 
review our model and the non-perturbative HEOM ap¬ 
proach. In Sec.HI, we present our accurate solutions re¬ 
lating to the PSB effect, those results include: HI.A. the 
physical picture of the PSB; HI.B. the gate voltage ma¬ 
nipulation of the PSB; and HI.C. the lift of the PSB by 
the dot-dot exchange interaction. In Sec.IV we give the 
summary of our work. 
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II. THEORY AND FORMULAS 


The HEOM is a general formula for the quantum open 
systems composed of three parts: the system (quantum 
dots here), the bath (two electrodes here) and the system- 
bath couplings. Let us introduce the total Hamiltonian 
(Anderson impurity model) for the DQD as follows, 

Ht = Hs + Hb + Hsb ( 1 ) 


where Hg is the Hamiltonian for the two coupled dots 


Hs=Y1 


U 


= 1.2c 

~\~t 'y ' (QijQ2o- + H.C.) (2) 


Ciaa^^ttia + 2 
i—l,2cr z—1.2cr 


here tia indicates the on-site energy of the electron with 
spin cr (cr =t,^^) on dot i{i = 1,2), and ciia cor¬ 
respond the creation and annihilation operators for an 
electron with spin tr. riia = is the electron num¬ 

ber operator of dot i, and U is the Coulomb interaction 
between electrons with spin a and a (opposite spin of tr) 
within one dot. t is the inter-dot coupling, determined by 
the overlapping integral of electron wave functions. H.c. 
stands for the Hermitian conjugate. 

For brevity, in what follows, we use the symbol /i to 
denote the electron orbital (including spin, space, etc.) 
in the system , i.e., fj, = The Hamiltonian of the 

electrodes is described as a noninteracting Fermi bath. 


from the basic path integral equations (influence func¬ 
tional theory) without take any approximations [Tsj l . The 
HFOM that governs the dynamics of the DQD takes the 
form of 










— I 


j3 ^3i---3n3 
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( 5 ) 


where the nth-order auxiliary density operator can 
be defined via auxiliary influence functional as 

( 6 ) 

with the reduced Liouville-space propagator, 

r'4’[t] f 

Jipo bo] 

5 [■(/>] is the classical action functional of the reduced sys¬ 
tem. The definition of the auxiliary influence functional 
together with its equations is referred to in [l5l| . 
We denote j = {ji ■••jn} and = 

{ji ■ • • jr-ijr-i-i • • • jn}, the action of superoperators 
respectively is 

APh +'y-) P'ij 


Hb — ^kad].^^dkfia (3) 

kp.a—L,R 

with Cka being the energy of an electron with wave vec¬ 
tor k in the a lead, and the d\^^{dkfia) corresponding 
creation (annihilation) operator for an electron with the 
a-reservoir state \k) of energy Cka- The dot-electrode 
coupling Hamiltonian is 

Hsb = ^[fl(t)df, + dlf^{t)] (4) 

in the bath interaction picture. Here, = 

is stochastic interactional 
operator and satisfies the Gauss statistics with takfi de¬ 
noting the transfer coupling matrix element. The influ¬ 
ence of electrodes on the dots is considered through the 
hybridization functions with a Lorentzian form, AQ,(a;) = 
'^J2k^akA*akti^i^-^ka) = AlT^/[2(a;-/ra)^-|-lT^], with 
A being the effective impurity-lead coupling strength, W 
being the band width, and p-a being the chemical poten¬ 
tials of the a lead. 

Obviously, Eq. © is a strongly correlated Hamiltonian 
with infinite degree of freedoms, which is hard to exactly 
solve by the the Schrodinger equation directly. Fortu¬ 
nately, we can derive the accurate HEOM for the re¬ 
duced density matrix (together with the auxiliary ones) 


(9) 

In this formalism, a° (a°) corresponds the creation (an¬ 
nihilation) operator for an electron with the p elec¬ 
tron orbital. The reduced system density operator 
p^^\t) = trs[ptotai(i)] and auxiliary density operators 
■ ■ ■ ) T} are the basic variables, here L 
denotes the terminal or truncated tier level. The Liou- 
villian of dots, £■ = ?i“^[i7sys, • ], contains the e-e in¬ 
teractions. The index j = {ppm) corresponds to the 
transfer of an electron to/from (o = +/—) the impu¬ 
rity state |(t), associated with the characteristic mem¬ 
ory time 7 “^. The correlation function C°^,^{t — t) = 
{faf,{t)fL{T))B follows immediately the time-reversal 
symmetry and detailed-balance relations. 

We set the initial total system at equilibrium where 
pa = = 0. The system will leave equilibrium after 

applying a voltage to the left (L) and right (R) leads, and 
there will be a current flowing into the a-lead Ia{t) 

Ia{t) = i'^As[pi^{t)d^ - dlp~^{t)] ( 10 ) 

Here, p\^^ = (p^^)^ is the first-tier auxiliary den¬ 

sity operator obtained by solving Eq.(5). Through 
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the extended Meier-Tannor parametrization method and 
multiple-frequency-dispersed hierarchy construction, we 
can achieved the closed HEOM formalism M- As a 
result, the current from L to R lead can be denoted 
I{t) = hit) = 

For evaluation of dynamical variables of the DQD sys¬ 
tem, we focus on correlation function between two ar¬ 
bitrary dynamical operators, CAB{t) = {A{t)B{0)) = 
trtotai[A(t)i?(0)p®^^^j(T)]. Here, the Heisenberg opera¬ 
tors and thermal equilibrium density operator 
are all defined in the total space. A linear re¬ 
sponse theory for quantum open systems [l^, 
has been established, based on which CAB{t) is re¬ 
trieved exactly within the HEOM framework. Let 
Cab{u}) = j J dte^^*CAB{t), which satisfies the de¬ 
tailed balance relation of Cba{—^) = 

The system spectral function is obtained as Jab{'-^) = 
^/dte*-*({i(t),H(0)}) = i (1 + 6-'^-/'=^^) Cab( w). 
With A = di_i and H = aj,, it recovers the spectral func¬ 
tion of the impurity state fi, i.e., A^{uj) = J. -t (w) = 
— ilmG^^(u;). Here, G^^(a;) is the retarded Green’s 
function. 

In our calculations, we treat the results as converging 
if the errors in numerical results of each element of the 
density matrix or the matrix of spectral function between 
the truncation L = N and L = -|- 1 are less than 5%, 

then sufficiently accurate current will be output. In the 
follow calculations we adopt L = 4 . 

The main advantages of the HEOM approach applying 
to the DQD systems are as follows: 1) the HEOM the¬ 
ory is established based on the Feynman-Vernon path- 
integral formalism, in which all the system-bath correla¬ 
tions are taken into consideration; 2) the HEOM method 
is nonperturbative. In principle, the HEOM formalism is 
formally exact for noninteracting electron reservoirs. It 
also resolve nonperturbatively the combined effects of e-e 
interactions; 3) the HEOM is a high-accuracy numerical 
approach. It has the ability to achieve the same level of 
accuracy as the latest high-level NRG method M- Its 
main disadvantage lies in the increasing computational 
cost as the system temperature decreases. 


III. RESULTS AND DISCUSSION 
A. Physical picture of Pauli spin blockade 

The spin degeneracy in Eq. ([IJ is not convenient for the 
detailed analysis of the PSB, thus we theoretically lift the 
spin degeneracy in dot I by means of a local magnetic 
field Bi applied onto it, with its direction paralleling to 
the down spins. Bi is chosen to be strong enough to make 
the energy level of up-spin electrons ei-;- much higher than 
eij,. In experiments, a local-like inhomogeneous Zeeman 
field can be archived by a novel split micromagnet [^ . 
By adjusting gate voltage Vi, we then set eij, at the equi¬ 
librium Fermi level Ep {Ep = pp = Pr at zero-bias). 


Under such condition, Ep coincides with the center of 
the peak of the transition for down-spin electrons in dot 
1 to jump from the zero occupied level to single occupied 
one, as the spectral function Ai 4 _(a;) shown in Fig.[T])a). 
This kind of single-occupied transition peak of Ai.|-(a;) is 
pushed far higher than Gb by the large Bi, which makes 
the up-spin current negligibly small, thus Ai.|-(a;) is not 
shown in figures [see Fig.|lja)-(c)] but its little contribu¬ 
tion to the total current is still counted [see Fig.[T][d)]. 

Without the coupling t between dot 1 and 2, the local 
Bi should take no direct effects on the spins in dot 2 . 
By adjusting gate voltage V 2 , we set the double-occupied 
transition (from the single occupied level to double occu¬ 
pied one) peaks of A 2 't'(a;) and ^ 24 ,( 0 ;) to coincide with 
Ep at t ~ 0, as shown in Fig.[T](a). Proceeding from 
this set-up, we gradually adjust the inter-dot coupling 
strength and other parameters to elaborate the physical 
picture of the PSB in details. 




A,ii«(ltenV) A^(1/nisVl 



FIG. 1 : (Color online). In the case of the spin non-degeneracy 
in dot 1, (a) equilibrium)!/ = 0) spectral function Ai^(a;), 
A2'[(ui) and ^24.(0;); (b) nonequilibrium spectral function 
Ai4(a;), A2^{ijj) and A2i{ijj) at positive bias(U = - 1-2 mV);(c) 
nonequilibrium spectral function Ai^(a;), A2-t-(w) and A24.(a;) 
at negative bias {V — —2 mV); and (d) the dependence of 
the occupation numbers Ni^, A24-, N21 as well as the total 
current I on the bias voltage V. The parameters adopted are 
Eit = £14. = £2t = £24, = —1.0, U\ = U 2 = 2.0, Vgi = — 2 . 5 , 
Vg 2 = 1.0, Wl = Wr = 4 . 0 , Fl = r_R = 2.0, AGsi = 1 . 5 , 
T = 0.1 and t = 0.001 (in unit of meV). 

We start our study on the PSB from the following pa¬ 
rameters: the potential energy £ 14 - = £ 14 ^ = £24 = £ 24 . = 
— 1.0 meV; on-site e-e interaction strength Ui = U 2 = 2.0 
meV; gate voltage Vgi = —2.5 meV and Vg 2 = 1.0 meV. 
The left and right electrodes are chosen to be symmetric, 
with their DOS being Lorentzian-type; bandwidth being 
Wl = Wr — 4.0 meV and the electrode-dot coupling 
being Tb = Tb = 2.0 meV. The Zeeman splitting energy 
caused by the local magnetic field Hi is AEri = 1.5 
meV; and the temperature T = 0.1 meV. The dot-dot 
coupling in Fig. [T] is near zero, t = 0.001 meV. 
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The spectral functions A{uj) shown in Fig.[lja) is cal¬ 
culated by HEOM at the equilibrium state (bias voltage 
E = 0), where only ^ 14 ^( 0 ;) is shown in dot 1. Fig.[T](b) 
and (c) show nonequilibrium spectral functions ^ 14 ,( 0 ;), 
A 2 ^{uj) and ^ 24 ,( 0 ;) at positive {V — +2 mV) and nega¬ 
tive {V = —2 mV) bias, respectively. Fig.[TDd) depicts 
the dependence of the occupation numbers N 2 -\, 

7 V 24 , as well as the total current I on the bias voltage 
V {N — V/1 — V curves), where N is obtained from 
the diagonal elements of the reduced density matrix, 
which is also equal to the weighted integral of corre¬ 
sponding A{uj) under the Fermi level at the steady state. 

As shown in Fig. [I] without the bias being applied 
{V = 0), Nil « 0.5, N 21 = N 21 « 0.7. The frac¬ 
tional charges results from the following two reasons: 
1) the value of the Fermi-Dirac function /(w) in the 
vicinity oi Ep continuously changes from 0 to 1 rather 
than a integer; and 2) the spectral function A{uj) near 
Ep shows a finite-width peak structure broadened by 
the dot-electrode interaction, thus the weighted inte¬ 
gral value below the Fermi surface is less than 1. The 
total electron number in the two-dot system is Qp = 
Vi|-|-Vi 4 ,-|-V 2 '|--l-V 2 ^ « 2.0 {Nil ^ 0 ), which remains al¬ 
most constant even in the non-equilibrium (steady-state) 
transport process if the external (electric and/or mag¬ 
netic) field is not too large. 

Then, the symmetrical positive {fip ^J■L + 

eV/2;nii Hr - eV/2) and negative {hl Hl - 
eV/2]HR Hr E eW^) bias are applied to the config¬ 
uration shown in Fig.[T](a). From the changes of A{uj), 
N and I with V shown in Fig.[T])b)-(d), we can see that 
the spectral function and occupation number of up-spin 
and down-spin electrons keep degenerate in dot 2 , and 
neither positive nor negative bias can lift the degeneracy. 
The reason lies in the very small coupling between two 
dots, t ^ 0 . 

The I—V curve [see Fig.[IJd)] shows distinct asymmet¬ 
ric behavior that the steady value of the current at V < 0 
is much larger than that at V > 0 , actually the former is 
almost twice of the latter. However, such asymmetry is 
not the same as the rectification since the positive steady 
current is not small enough (comparing to the negative 
one ) to define a blockage effect. That result can be ex¬ 
plained as follows: 

In principle, the dot-to-dot electron transfer can induce 
an antiferromagnetic exchange between them with the 
strength Jap ^ At'^/U. Enough strong Jap will lock the 
ground state of the isolate double dot system into a spin 
singlet state S'(l, 1) whose energy is lower than the triplet 
state T(l, 1 ) in the order of Jaf- In Fig.[Tl t ^ 0 means 
Jaf ^ 0, thus states of S'(l,l) and T(l,l) are nearly 
degenerate at V = 0. 

When positive bias applied, hr < ns shown in 
Fig.HTd), N 21 and N 21 decreases gradually with hr de¬ 
creasing . The relation of N 21 ~ V 24 , at V = 0 will keep 
unchanged at V > 0 , and the number will reach a steady 
value after V > 1 mV, N 21 ~ ~ 0.5. At the same 


time, Nil increases gradually with hl{hl > 0 ) increas¬ 
ing and tends to 1.0 after V > 1 mV. The nonequilibrium 
spectral function at V = +2 mV in Fig.jTJb) further con¬ 
firms the above process. It means that the states S'(l, 1) 
and 'r(l,l) remain degenerate at V > 0 , and each of 
them has 50% probability after V > 1 mV. It is believed 
that electron spin is conserved through direct hopping 
process,thus the initial state T(l, 1) can only transfer to 
T(0, 2), if one electron is driven from dot 1 to dot 2 by the 
positive bias. However, the state r(0,2) is not permit¬ 
ted to exist due to the Pauli’s exclusion principle (the 
exception from the excess freedom such as the orbital 
or valley not considered in the present work). Therefore, 
only half of the total initial states [S'(l, 1 )] can contribute 
to the transport current at V > 0 via the transition 
^(1,1) ^5(0,2). 

As shown in Fig.[T][d), the negative bias makes hr > 0i 
then ^ 24 - and N 21 increases gradually with hr increas¬ 
ing but remains ^ 24 - ~ A^ 24 , which is close to 1.0 after 
V < — 1 mV. Meanwhile, Nii decreases gradually with 
Hl{hl < 0) decreasing and approaches to 0 after V < — 1 
mV. The nonequilibrium spectral function at V = —2 
mV in Fig.[T](c) further confirms the above process. It 
suggests that the negative bias will stabilize the S'(0,2) 
state which 100 % contributes to the current, since the 
Pauli’s exclusion principle takes no effect on the transi¬ 
tion from ^(O, 2) to S'(l,l). That argument is verified 
by the I — V curve in Fig.[ljd), which shows that the 
steady current at V < 0 is almost twice larger than that 
at V > 0 . 



I ■ I ■ 1 ■ I ■ I 

- 2-1012 
Bias V(mV) 


FIG. 2 : (Color online). In the case of the spin non-degeneracy 
in dot 1, (a) equilibrium)!/ = 0) spectral function ^14(0;), 
A24-(a;) and A2i{ij); (b) nonequilibrium spectral function 
Ai4(a;), A2-[{ijj) and A2i{ijj) at positive bias)!/ = - 1-2 mV);(c) 
nonequilibrium spectral function ^14(0;), A24-(a;) and ^24(0;) 
at negative bias {V — —2 mV); and (d) the dependence of 
the occupation numbers V14, ^24-, A24 as well as the total 
current I on the bias voltage V. The parameters adopted are 
the same as those in Fig.[Tl except t — 0.08 meV here. 

We then increase the inter-dot coupling strength t from 
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0.001 to 0.08 meV and keep other parameters unchanged 
for the purpose of comparison. The calculated A(w, V = 
0), A(uj, V = +2 meV), A(uj, V = —2 meV) and the N — 
V/I — V curve are depicted in Fig.[2Ua)-(d) respectively. 
From Fig.[2ja) for the case of F = 0, we can see that the 
degeneracy of S'(l, 1) and '^(l, 1) has been lifted by the 
t-induced antiferromagnetic interaction Jaf{'^ At^/U). 
Since the energy of <5'(1,1) is lower than that of T(l, 1), 
A 2 -f{uj) moves downward and N 2 -\- increases from ^ 0.7 
at t ~ 0 to ~ 0.72. Meanwhile, N 21 decreases from ^ 0.7 
to ^ 0.68, as shown in Fig.[5][d). 

The PSB will take place when the positive bias is ap¬ 
plied to the set-up shown in Fig.[^a). At F >0, fV2't' 
and N 21 changes to different directions [see Fig.[5][d)] in¬ 
stead of synchronous varying at t ^ 0 [cf. Fig.[T][d)]. As 
shown in Fig.[2][d), A^2t decreases rapidly as F positively 
increasing and tends to 0 after F > 1.0 mV; whereas 
N 21 increases and approaches 1.0 after F > l.OF mV. 
As for Nil, if gradually increases from 0.5 to 1.0 with 
increasing F, much like the change of N 21 . As a result of 
above changes, only state T(l, 1) is retained under pos¬ 
itive bias and 5'(1,1) will no longer exist after F > 1.0 
mV. As already explained, T(l,l) can not transfer to 
T(0, 2) to create any current due to the Pauli’s exclu¬ 
sion principle, thus PSB occurs naturally. The forbid¬ 
den electron transition from dot 1 to 2 outputs near zero 
current, as the I — V curve shows in Fig.[21[d). That 
is exactly the PSB effect observed in experiments at a 
moderate coupling strength t, and the finite current in 
the interval of F S [0,1.0 mV] is the so-called leakage 
current corresponding to the process of S'(l, 1) being de¬ 
pleted gradually. After F > 1.0 mV, the current enters 
its total-blocked zone with a near zero value. 

When negative bias applied as shown in Fig.[5][d), with 
the increase of fiR, N 2 -\- increases gradually and becomes 
saturated at iV2-|- ~ 1.0 after the bias F < —1.0 mV, while 
N 21 decreases and keeps at about 0.5 + 5 (5 ~ 0.13) after 
F < —1.0 mV. At the same time, Ni^ decreases with 
decreasing and maintains Ni^ Ki 0.5 — 5 after F < —1.0 
mV. Although Ni^ > 0 in this case, it won’t cause any 
PSB effect of the down-spin electrons, for the reason that 
the below the Fermi surface can provide enough space 
to accept the electrons transferring from dot 2 to 1, as 
the nonequilibrium spectral function shown in Fig. He). 
As a consequence, the considerable current is output in 
the I — V curve. It should be noted that the unit of 
current in Fig.[^d) is pA instead of nA in Fig.[ljd). 

Summarizing Fig.[I]and[21 our theory appropriately de¬ 
scribes the physical mechanism and picture of the PSB 
effect in DQD systems, by means of a strong local mag¬ 
netic field applied onto dot 1 to lift its spin degener¬ 
acy. We are now on the position to elucidate the PSB 
effect under more general conditions. In Fig. [31 we depict 
the N — V/I — V curves at t — 0.001 meV [Fig.j^a)] 
and 0.08 meV [Fig.|3l[b)], where the other parameters 
are the same as those in Fig.jT] and Fig. [3] accordingly, 
except that the local magnetic field is absent in both 
cases, i.e. Bi = 0. By comparing the I — V curve in 


Fig.|3][a) to that in Fig.|T][d), we can see that the ad¬ 
ditional up-spin channel barely changes the current at 
F > 0 but will increase that at F < 0 if the dot-dot 
coupling is very weak. The reason for the former lies 
in the equally dividing of Ni^ in Fig.[T][d) by Ni^ and 
Nil Fig-Hs-) the current changing little. The rea¬ 
son for the latter is that the up- and down-spin channels 
contribute to the current independently in the limit of 
t ~ 0, thus the current will be enhanced by additional 
channels. If distinct PSB effect occurs, the situation will 
become much different. As indicated by the comparison 
of the I — V curve in Fig.|3][b) to that in Fig.|2][d), the 
additional up-spin channel hardly changes the current 
either in the PSB region (F > 0) or in the conductive 
region (F < 0). By analyzing the corresponding N — V 
curves and nonequilibrium spectral functions (the figures 
not shown), we find in the PSB region, the probability 
of r(4., 4 .) in Fig.|3][d) is equally divided into r(t, t) and 
T(4,, 4 .) in Fig.inDb), and the current keeps its value. In 
the conductive region, the single-spin transport channel 
5'(0,t4-) S'(4., t) is equally divided by the degenerate 

double-spin ones, S'(0,t4') ^ ^‘S'(4'i t) + i); and the 

total current remain unchanged. 




FIG. 3: (Color online). Without the local magnetic 
field applied {Bi — 0), the dependence of the occupation 
numbersAi-I", ^14., ^24-, and N 21 as well as the total current I 
on the bias voltage F at t = 0.001 meV (a) and 0.08 meV (b). 
The other parameters in (a) and (b) are the same as those in 
Fig. [T] and Fig.[2] accordingly. 


B. Gate voltage modulation on PSB 

In experiments on PSB, the gate voltages Fi and F2 are 
two important parameters which respectively manipulate 
the on-site energy of dot 1 and dot 2. We thus theoreti- 
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cally investigate the variation of the PSB with them and 
summarize the results in Fig. 01 where the dot-dot cou¬ 
pling strength is chosen as t = 0.05 meV, a small value 
but large enough to induce the PSB effect, with the pur¬ 
pose of making the boundaries shown in Fig. 0] clear and 
distinguishable. Fig.0] (a) depicts the positive current at 
V = OA mV changing with respect to the parameters in 
the Vi — V 2 plane, in the form of the 3D colormap surface 
image together with the 2D bottom contour projection. 
The other parameters are the same as those in Fig.[3l ex¬ 
cept both Vi and V 2 are variables now within the range 
of [—U, U]{U = 2 meV). In the 2D projection image, we 
schematically mark off the boundary of the stability dia¬ 
grams by dotted lines. Actually, the quadrangles shown 
in the figure should changes to hexagon at finite t, how¬ 
ever, the boundary line is hard to accurately determine in 
theory. The schematic stability diagrams shown in Fig.0] 
is just for reference purposes. 



FIG. 4 : (Color online). 3 D colormap surface image with the 
2 D bottom contour projection of the current (absolute value) 
as a function of the gate voltage Vi and V2. The bias voltages 
used are (a) V = 0.4 mV, and (b) V = — 0.4 mV. In the 2 D 
projection image, the dotted lines schematically mark off the 
boundary of the stability diagrams. 

Fig.0] clearly shows how (14,14) modulate (^ 1 ,^ 2 ) 
and the PSB effect. The DQD system remains the stabil¬ 
ity charge-occupied state (1,1) within —U/2 <Vi< U12 
and —U/2 < V 2 < D/2, and no current occurs in the 
center of this area due to the Coulomb blockade. How¬ 
ever, finite current may be output at the four top corners 
of ( 1 , 1 ) state, resulting from the charge transferring in 
and out, which can be listed as ( 0 , 0 ) ( 1 , 1 ), ( 2 , 0 ) ^ 

( 1 , 1 ), ( 1 , 1 ) ( 2 , 2 ), and ( 1 , 1 ) ^ ( 0 , 2 ) (clockwise from 


lower left). By referring the figure, one can see that the 
current at the corners of ( 0 , 0 ) ^ (1, 1 ) and ( 1 ,1) ( 2 , 2 ) 

is symmetric about the bias voltage with no PSB effect 
occurring. In addition, the current at above two corners 
also shows symmetric behaviors along the diagonal line 
(14 = 14 ), which comes from the electron-hole symmetry 
satisfied by our Hamiltonian, Eq. ©■ As shown in Fig. 01 
current at the corners of ( 2 , 0 ) ^ ( 1 , 1 ) and ( 1 , 1 ) ( 0 , 2 ) 

exhibits rectifying characters about positive and negative 
bias. The former case has been elaborated in Fig.|T]to01 
and the small current (|/| < 1 pA) shown in Fig. 0Kb) 
corresponds to the leakage current at low bias shown in 
Fig. [3] The latter case of (1,1) ^ (0,2) is very similar 
except that the PSB effect takes place at the positive 
bias. 

By referring Fig. 01 one can see that large current 
more than 4 pA can occur at the top corners of 
(1,1) state, namely the center points of transition be¬ 
tween (1,1) and other stability states. In Fig.0Ka), 
those points correspond to (—D/2,—D/2), (D/2,—D/2) 
and {U/2, U/2) in the (14,14) parameter plane, versus 
{-U/2, -U/2), {-U/2, U/2) and (D/2, D/2) in Fig.0Kb). 
What special about those points is that one quantum 
transition peak in dot 1 will resonate with another one 
in dot 2 at the Fermi surface, as the nonequilibrium spec¬ 
tral functions shown in Fig. 01 where Fig.[5Ka) to (d) cor¬ 
responds to the point (—D/2,—D/2) to (—D/2, D/2) in 
the clockwise order. Taking Fig.jSKa) as an example, we 
can see that the resonance on point (—D/2, —D/2) takes 
place between the transition of (0 —> 1) in dot 1 and that 
of (0 —1) in dot 2 at the Fermi surface, which induces a 
very large current, / ^ 4.2 pA. If (14,14) deviates those 
points parallel to the 14 = 14 diagonal line, the resonance 
between the transition peaks will still exist, but no longer 
coincide with the Fermi surface. As a consequence, the 
current will gradually decrease into the Coulomb block¬ 
ade region, after some peak-like structures, as shown in 
the 3D colormap surface image in Fig.01 If (14,14) devi¬ 
ates the four top corners of (1,1) state along any direction 
of 14 yf 14 , none of the resonance will survive, and then 
the current will decays to very small value {I < 0.5 pA) 
quickly. 

Although both peaks of the transition of (0 —?> 1) in 
dot 1 and that of (1 —>- 2 ) in dot 2 are coincide with the 
Fermi surface on the point (—D/2, D/2), the PSB pro¬ 
hibits the resonance between them as shown in Fig.jSKd), 
and only very small leakage current (/ < 1 pA) occurs 
in the vicinity of (—D/2, D/2) [see Fig.0Ka)]. When the 
direction of the bias is reversed from V > 0 to V <0, 
{—U/2, U/2) will change from a PSB point to a resonance 
one, and then very large current occurs at this point. Ac¬ 
cordingly, the PSB point moves to (D/2, — D/2), as shown 
in Fig.HKb). 

As indicated in Fig.01 our theoretical results of the 
gate voltage modulation on PSB is qualitatively consis¬ 
tent with the experimental measurements. The shape 
of conduct regions is also similar to the triangles ob¬ 
served in experiments, but not exactly the same. The 
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one Jaf induced by t, that is Jt ~ + 4t^/C/. In or¬ 

der to highlight the role of J, we thus choose a relatively 
small t {t = 0.05 meV) at fixed U = 2.0 meV to produce 
a very small Jaf {Jaf ^ 0.005 meV). 

Fig.[S] shows our results of gate voltage modulation 
(bias V = 0.4 mV) on PSB at various J with dif¬ 
ferent signs and values, where Fig.[5Da), (b), (c) and 
(d) correspond to J = —0.32, —0.08, 0.08 and 0.32 
meV, respectively. By referring Fig.jSJ one can see 
that J substantially affects the current near the points 
(Vi, V 2 ) = {-U/2, U/2) and (C//2, -17/2) where the PSB 
effects take place, while it slightly dose to the current 
near {—U/2, —U/2) and {U/2, U/2) where no PSB effects. 
It indicates that the exchange interaction can directly 
change the characters of the PSB. 


FIG. 5: (Color online). Non-equilibrium spectral functions at 
four top corners of (1,1) state in Fig.|4] (a), (b), (c) and (d) 
respectively represent the point {—U/2, —U/2), {U/2, —U/2), 
{U/2, U/2) and {—U/2, U/2) in the {Vi, V 2 ) parameter plane. 


difference may come from the reason that the Anderson 
multi-impurity model can not describe all the details in 
experiments. For instance, when the gate voltage Vi on 
dot 1 changes, it has been confirmed by experiments that 
the effects on dot 2 are induced not only through the di¬ 
rect coupling t but also through a capacitive coupling. 
We believe the former has been well described in our 
theory, but the latter has not yet. 
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C. Lift of PSB by exchange interaction 

Now we extend the Anderson two-impurity model to 
adding the term of dot-dot exchange interaction, which 
describes the coupling between the local spins of two dots, 
with the Hamiltonian as follows: 

He = JSi ■ S 2 ( 11 ) 

where Si is the spin operator of dot i{i = 1,2). J 
is coupling strength between spins in different dots, 
which could be positive(antiferromagnetic) or nega- 
tive(ferromagnetic). 

Fundamentally, J originates from the exchange term of 
the e-e interaction (potential Energy) between two dots, 
thus plays an equal important role as the kinetic energy 
t. The manipulation of J is a fascinating issue closely rel¬ 
evant to quantum information . Despite the practical 
difficulties, experiments may achieve this goal indirectly, 
for example, N. J. Craig et al. have demonstrated the 
control of the strength and the sign of J between two 
dots coupled through an open conducting region [ 2 ^. In 
the present work, we investigate the effect of J on PSB. 
When t 0, the total exchange interaction Jt princi¬ 
pally equals to the sum of J and the antiferromagnetic 


FIG. 6: (Color online). 2D contour plot of the current (bias 
V=0.4 mV) as a function of the gate voltage Vi and V 2 at 
various J. (a), (b), (c) and (d) correspond to J = —0.32, 
—0.08, 0.08 and 0.32 meV, respectively. 

Let us focus on the change of the leakage current of 
the PSB around {—U/2,U/2). It corresponds to a very 
small value at J = —0.32 meV [Fig.jSKa)], which indicates 
that the ferromagnetic exchange interaction tends to en¬ 
hance the PSB effect. Generally speaking, the leakage 
current will increase with the increasing of the algebra 
value of J, as shown in Fig. [51 For example, the leak 
current increases from near zero to ~ 0.5 nA [Fig. [6Kb)] 
as J increasing from —0.32 to —0.08 meV. Changing 
the sign of J still maintains this kind of tendency, and 
the antiferromagnetic exchange interaction seems to sup¬ 
press the PSB. When J continually increases to 0.08 meV 
[Fig.lHKc)], the leakage current is clearly visible and in the 
range of I ~ 2 nA. If J positively increases to a enough 
large value, e.g. J = 0.32 meV [Fig.lHKd)], the leakage 
current will increase distinctly, even to the same order of 
magnitude (~ 4 nA) as the conductive current. In this 
case, we argue that the PSB has been lifted by strong 
antiferromagnetic exchange interaction. 

The lift of the PSB is a significant feature of Fig. [51 


























































which deserves careful study to reveal its mechanism. We 
thus apply the local magnetic field Bi again onto dot 1 
to lift is spin degeneracy and investigate the transport 
of down-spin electrons in the PSB region. Fig. [7] shows 
the nonequilibrium spectral functions at the point of 
(Vi,V 2 ) = (—17/2,17/2), with the Zeeman energy caused 
by Bi being AEbi = 1-5 meV and other parameters be¬ 
ing the same as those in Fig. [01 Fig.|7Ka)-(d) corresponds 
to Fig.| 6 Ka)-(d), respectively. 

At J = —0.32 meV, the dot-dot exchange interaction 
is ferromagnetic, which makes the energy level of T(l, 1) 
lower than 5(1, 1 ). As a consequence, the single occupied 
transition peak of down-spin in dot 2 is lower than that 
of up-spin, since that peak in dot 1 has been locked as 
down-spin one at the Fermi surface. It means that the 
down-spin electron in dot 2 has almost fully occupied the 
single level under the Fermi surface, which will prevent 
the hopping of electrons with the same spin and thus 
enhance the PSB effects. Continuously increasing the al¬ 
gebra value of J to —0.08 meV will not essentially change 
above process, and only slightly increase the weight of the 
down-spin holes in the double occupied transition peak 
at the Fermi surface, which makes the leakage current 
increase to a small nonzero value, as shown in Fig.jOKb). 



holes in the double occupied transition peak increases fur¬ 
ther, offering more space for down-spin electrons to trans¬ 
fer from dot 1 , and consequently increases the leak cur¬ 
rent value. Since J is small in Fig.ITJc), the leak current 
can only increase to a relatively small value, as shown in 
Fig.jSKc). In Fig.jTKd), large enough J (J = 0.32 meV) 
makes the weight of down-spin hole increase remarkably, 
thus induces the large leak current shown in Fig.jSKd) 
which lifts the PSB effect. By referring Fig.|6]and|71 we 
can see that the lift of PSB by J is a continuous process 
instead of a sudden change. 





FIG. 7: (Color online). In the case of the spin non-degeneracy 
in dot 1, nonequilibrium (bias V=0.4 mV) spectral function 
Ai^(a;), A 2 t(a;) and ^ 2 ^( 0 ;) at (^ 1 ,^ 2 ) = (-[//2,17/2). (a), 
(b), (c) and (d) correspond to J = —0.32, —0.08, 0.08 and 
0.32 meV, respectively. 

When J > 0, the spin exchange interaction is anti¬ 
ferromagnetic, which makes the energy level of 5(1,1) 
is lower than that of T(l, 1). This case should help elec¬ 
trons in dot 1 transfer to dot 2 by means of the transition 
5(1,1) —>■ 5(0, 2), so the leak current will increase, which 
has been confirmed by the change of the spectral func¬ 
tions shown in Fig.ITJc) and Fig.jTKd). We can see that 
after the sign change of J, the weight of the down-spin 


FIG. 8: (Color online). In the case of the spin non-degeneracy 
in dot 1, the dependence of the occupation numbers Ni^, N 2 ^, 
and N 2 i as well as the total current I on the bias voltage V 
at J = —0.32 meV (a) and 0.32 meV (b). 

To further verify the mechanism of the lift of PSB by 
the exchange interaction, we calculate the N — Vj I — V 
curves at different J and summarize the results in Fig. [51 
where Fig.jSKa) corresponds to J = —0.32 meV and 
Fig.jHKb) to J = 0.32 meV. By comparing those two fig¬ 
ures in Fig.jSl one can see that A^ij, (J = —0.32 meV)~ 
-^ 14 , (<7 = 0.32 meV) at bias V = 0, but both N 21 and N 2 ^ 
are distinctly different at J = ±0.32 meV. Specifically, 
when J = —0.32 meV and V = 0, N 21 > N 2 f, which 
means that the single occupation of down-spin electron 
in dot 2 is sufficient and the weight of hole with the same 
spin in double occupation is small, thus the transfer of 
down-spin electrons from dot 1 to 2 will be blocked. With 
the bias positively increasing, as shown in Fig.jSKa ), 
gradually increases and approaches to 1.0 after F > 1.0 
mV, as a consequence, the PSB effect is enhanced and 
the current decreases to a near zero value. Negatively 
increasing bias will produce a steady current at about 
—4.0 nA after F < — 1 mV as shown in the same figure. 
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On the other hand, when J = 0.32 meV and V = 0, 
^21 < N 2 -f, which means that the single occupation of 
down-spin electron in dot 2 is small and the weight of 
hole with the same spin in double occupation is large, 
thus the PSB of down-spin electrons from dot 1 to 2 will 
be lifted. With the bias positively increasing, as shown 
in Fig.IHUb), 7 V 24 , slowly increases and then stabilises at 
0.9 after V > 1.0 mV. It suggests that the weight of hole 
in double occupation keeps finite, which will induce con¬ 
siderable leakage current even at large positive V. Neg¬ 
atively increasing bias will produce a steady current at 
about —6.0 nA after V < — 1 mV, much larger than that 
at J = —0.32 meV [cf. Fig.|8] (a) and (b)]. The rea¬ 
son may lie in the fact that antiferromagnetic interaction 
tends to form S'(l, 1) which is in favor of the transition 
of 5(0,2) ^ 5(1,1). 

We are now on the position to compare the different 
roles of J and t on PSB. For this purpose, we define a 
physical quantity Iq '■ II as the ratio between the con¬ 
ductive current /c at V = —2.0 mV and the leakage one 
/l at V = 2.0 mV at the point of (Fi, V 2 ) = (—17/2,17/2), 
and then summarize the dependence of /q : /l on J in 
Fig. [Altogether with Iq : II on t in its insert. In principle, 
the larger the ratio Iq '■ Il is, the smaller the leakage cur¬ 
rent, and the more distinct the PSB effect, in case that 
Ic does not change significantly. As shown in the figure, 
the ratio Ic '■ II can reach 21 at J = —0.32 meV, indicat¬ 
ing a well-defined PSB effect, as shown in Fig. IlKa) and 
Fig.lSKa). Positively increasing J will induce a contin¬ 
uously decreasing of /c : II, smoothly passing through 
the zero point, /c:7i;,~15atJ = 0 , and finally ap¬ 
proaching a small value Ic '■ Il i si, J = 0.32 meV, 
indicating a total lift of the PSB effect. At J < —0.2 
meV and J > 0.2 meV, one can see that the decrease 
oi Ic '■ IL with J exhibits a near-linear behavior, but a 
nonlinear swell around the zero point is clearly shown in 
the figure. It suggests that J = 0 is a crossover point 
from the appearance to the lift of the PSB effect (at a 
fixed small t). 



J (meV) 


FIG. 9: (Color online). The dependence of the ratio Ic{V = 
—2.0mV) : 1l{V = 2.0mV) on J at the point of (Vi,V 2 ) = 
(—1//2, [7/2). The insert shows the dependence of the ratio 
Ic ■ II on t. 


By referring the insert of Fig.|9l one can see that the 
dependence of Ic : II on t exhibits a peak structure. 
More specifically, starting from the non-PSB point at t ^ 
0 [see Fig. |T] and Fig.jAKa)], increasing t will firstly result 
in a rapid increase of /q : Jl to a value about 16 at 
t = 0.8, denoting a well-defined PSB effect, as shown in 
Fig. [Hand Fig.|3Kb). By further increasing t {t > 0.8), the 
Ic ■ II will not increase any more, but slowly decrease to 
a small value about 5 at t = 0.3. It represents that large 
t does not mean distinct PSB effects which should only 
takes place at moderate t. That result is consistent with 
the experimental measurements in literatures 1 ^. 


IV. SUMMARY 

In summarize, we systematically investigate the Pauli 
spin blockade in double quantum dot systems. We start 
from the Anderson multiple impurity model to describe 
the system, fully considering the electron-electron inter¬ 
action and the dot-electrode couplings. By using the hier¬ 
archical equations of motion approach, we deal with this 
quantum model non-perturbatively to accurately obtain 
the equilibrium and nonequilibrium spectral functions, 
occupation numbers and current, etc. 

By means of a strong local magnetic field applied onto 
dot 1 to lift its spin degeneracy, our theory appropriately 
describes the physical mechanism and picture of the Pauli 
spin blockade effect in double quantum dot systems, fol¬ 
lowed by a general discussion without the local field ap¬ 
plied. Then, the gate voltage manipulation of the spin 
blockade is elaborated in detail by our theory. Our re¬ 
sults are proved to be qualitatively consistent with the 
experimental measurements. 

We further extend the Anderson multiple impurity 
model to involve the dot-dot exchange coupling, and care¬ 
fully study its effect on the spin blockade by changing the 
strength and the sign of the coupling. It is found that the 
ferromagnetic exchange interaction tends to enhance the 
spin blockade, while the antiferromagnetic one to sup¬ 
press it. What is more, the Pauli spin blockade effect 
may be lifted by the strong antiferromagnetic exchange 
coupling. 
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